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Abstract The measurement method of mode II fracture 
toughness-crack propagation length relation (i.e., the resis¬ 
tance curve, or R-curve) was examined by end-notched 
flexure tests on sitka spruce (Picea sitchensis Carr.). The 
tests were conducted by varying the span/depth ratios under 
the constant loading point displacement condition. The 
fracture toughness was measured from the load-crack shear 
displacement (CSD) and load-longitudinal strain relations. 
The crack length was determined by a combination of load- 
CSD and load-strain compliances and Williams’s end 
correction theory, as well as the observation of crack propa¬ 
gation. When the specimen had an appropriate span/depth 
ratio, the fracture toughness and crack propagation length 
were measured from the load-CSD compliance and com¬ 
bined load-CSD and load-strain compliances, respectively, 
and the R-curve could be determined properly under the 
constant loading point displacement condition. 

Key words Mode II fracture toughness ■ Crack length • R- 
curve • End-notched flexure test 


Introduction 

Wooden beams usually contain natural fractures whose 
planes parallel the longitudinal direction; and bending load 
is often applied to beams with these fractures. With this 
loading pattern, fracture often propagates by in-plane shear 
mode (mode II) or a combination of crack opening (mode 
I) and in-plane shear modes. Therefore, it is important to 
measure the mode II fracture toughness properly. When 
measuring the mode II fracture toughness of advanced com¬ 
posites such as carbon fiber-reinforced plastics (CFRP), the 
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end-notched flexure (ENF) test is commonly performed; its 
techniques were developed over two decades. 15 

In previous studies the ENF test method was examined 
using western hemlock, and it was found that the testing 
method could be applied to measuring the fracture tough¬ 
ness of wood under certain testing conditions. 6-7 The ENF 
tests were conducted under the loading point displacement 
control condition. Although the crack propagation is not 
always stable with this loading condition, a servo-controlled 
testing machine is not required and plural fracture 
toughnesses can be obtained when the initial crack length is 
l/ 3 v/3 (= 0.693) times longer than the half-span. Therefore, 
the fracture toughness-crack propagation length relation 
(i.e., the resistance curve, or R-curve) can be obtained even 
in this loading condition. As for the mode II fracture tough¬ 
ness of wood, however, there were few studies of fracture 
toughness by the R-curve because of the difficulty of evalu¬ 
ating the crack length properly. With the ENF test, the 
crack often propagates obscurely; and the crack length, 
which has a serious influence on the R-curve, is difficult to 
measure by observation. Kageyama and colleagues pro¬ 
posed the determination method for crack length without 
observing the crack propagation based on “Williams's end 
correction theory.” 3 ' 4 Nevertheless, Williams’s theory re¬ 
quires the elastic constants of the specimen, which should 
be determined independently of the ENF tests; and there is 
a concern that the variation in elastic constants inevitably 
contained in wood cause an incorrect evaluation of crack 
length. Here a method for determining the mode II R-curve 
of wood is proposed, and the validity of this method was 
examined. 


Theories 

Derivation of energy release rate by ENF test 

As in Fig. 1, an ENF specimen with a crack length of a is 
loaded by R at the center of the span of 2 L. The equation 
for bending is given by the elementary beam theory as 
follows. 
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(0 < x < a) 


(a < x < L) 
(L < jc < 2 L) 


( 1 ) 


where E x is Young’s modulus in the longitudinal axis, and / 
is the second moment of cross-sectional area in the crack- 
free region. Solving this equation, we can obtain the deflec¬ 
tion at the loading point v as: 


P(2L 3 + 3a 3 ) 

12^7 


( 2 ) 


The loading point compliance C L is defined as follows. 


_ v _ 2 L 3 + 3o 3 

L ~~ P 12 E X I ^ 

Hence, the energy release rate in the in-plane shear mode 
(mode II) , G„, is given as 1 ' 3 ’ 4 : 


r = — = 9C l pV 

11 2b da 2b{2L 3 + 3a 3 ) ^ 

where b is the width of the beam. When dGJda is nega¬ 
tive, the crack propagates stably. For stable crack pro¬ 
pagation, the crack length a should be larger than L/ 3 V3 
under the constant loading point displacement condition. 1 

Kageyama and colleagues proposed a method for deriv¬ 
ing the energy release rate by measuring the crack shear 
displacement (CSD) instead of measuring the loading point 
displacement as in Fig. I. 3 The CSD, denoted as <5, is ob¬ 
tained by integrating the strain along the longitudinal axis 
from the lower end corner of the upper beam to the crack 
tip: it is given as follows. 


d = 


ha 2 _ ha 2 
8 E X I' ~ ~EJ 


(5) 


where h is the upper beam height of cracked region, and I' 
is the second moment of cross-sectional area in the cracked 


region (which equals 118). The CSD compliance C c is writ¬ 
ten as follows. 



( 6 ) 


From Eqs. (3)-(6), the energy release rate can be formu¬ 
lated as 3 ' 5 : 


_ P 2 dC L _ P 2 2a- _ 3 P 2 C c 
11 ~ 2b da ~ 2b' 4E X I ~ 8 bh ^ 

There is a concern that the fracture toughness obtained 
from Eq. (4) or (7) is influenced by the additional deflection 
caused by the shearing force that implicitly occurred in the 
beam. 7 When the longitudinal strain at the most outer plane 
is measured (Fig. 1), however, the influence of shearing 
force might be effectively reduced. As in Fig. 1, the longitu¬ 
dinal strain at the midspan e x is represented as follows. 


Mh _ PLh 

~EJ ~ 2 EJ 


( 8 ) 


where M is the bending moment at the midspan. When the 
load-strain compliance is defined as C s , it is represented as: 



Lh 

2EJ 


(9) 


From Eqs. (3), (4), and (9), the energy release rate is given 
as: 


P 2 dC L _ 3 P 2 a 2 C s 
2b da 4 Lbh 


( 10 ) 


Evaluation of crack length 


In the ENF test the crack propagates without a clear open¬ 
ing, and the difficulty evaluating the crack length causes the 
impropriety of the A'-curve even when the fracture tough¬ 
ness is properly measured. To reduce this difficulty, the 
crack length must be determined without observation. One 
method for evaluating the crack length is based on 
Williams's end correction theory. According to that theory, 
the crack length is regarded as being -/h longer than the real 
length, so the load-CSD compliance C c is represented as 
follows. 2 


c = h(a + xhf 
C EJ 

where % is the function of elastic constants as: 


1 = 



3-2 


4 E *E y 


v 4.8 G xy 


^E 


y J 


( 11 ) 


( 12 ) 


where E y is Young’s modulus in the direction of thickness, 
and G xy is the shear modulus. The crack length evaluated by 
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the theory is defined as a w and is determined from Eq. (11) 


= /; 


~X + 


2 bE r C r 


(13) 


With Williams’s end correction theory, which was origi¬ 
nally applied to measuring the fracture toughness of CFRP 
by Kageyama and colleagues, 3 measurements of elastic con¬ 
stants are required independently of the ENF tests. Never¬ 
theless, the elastic constants of wood show larger variations 
than those of CFRP, and there is a concern that these varia¬ 
tions cause the impropriety of R-curve. Additionally, it is 
inconvenient to measure the elastic constants indepen¬ 
dently of the ENF tests. By combining the compliances, 
however, the crack length during propagation can be calcu¬ 
lated without measuring the elastic constants. From Eqs. 
(3), (6), and (9), the compliance ratios r LS , r cs , and r CL are 
defined by C L /C S , C c /C s , and C c /C L , respectively, as follows. 



r _ y 

L 2 

2 + 3 - 

C L _ 

L uj j 

Cs 

6 h 



Normalized crack length a/L 


Fig. 2. Variances of normalized compliance ratios by the crack 
propagation 


Experiment 

Materials 


r cs ~ 



Cc 

C L 



L 



(15) 


(16) 


The ratios at a — 0.7 L, which is close to the lower boundary 
for stable crack propagation, are defined as r° LS , r° cs , and r° ci _. 
Figure 2 shows the variances of normalized compliance ra¬ 
tios r LS /r^, r cs lr° cs , and r CL /r° CL by crack propagation. As 
shown in Fig. 2, the variance of r cs /r° cs is larger than that of 
the others, which indicates that the value of r cs is more 
sensitive to the change in crack length than the others. 
When evaluating the crack propagation length, therefore, it 
would be effective to use the combination of load-CSD and 
load-strain compliances and the crack length estimated by 
the combination of compliances, which is defined as a c , 
derived from Eq. (15) as: 




(17) 


In this experiment, the crack length was evaluated by a, a w , 
and a c . The R-curves obtained by the various crack length 
evaluations were compared with each other. 


Sitka spruce (Picea sitchensis Carr.), with a density of 0.42g/ 
cm 3 , was used for the specimens. Specimens were condi¬ 
tioned at 20°C and 65% relative humidity (RH) before and 
during the tests. Five specimens were used for each testing 
condition. 


ENF tests 

A beam specimen was cut with the dimensions of 15 X 
15 mm in the radial and tangential directions, respectively, 
and 250-550 mm at intervals of 50 mm in the longitudinal 
direction. The span length 2 L varied as 150-450 mm, which 
was 100 mm shorter than the total length. The crack was 
initially cut in the longitudinal-tangential plane by a band 
saw with a thickness of 1mm; it was then extended by a 
razor blade. The initial crack length a Q was 0.7 times the 
half-span L. A strain-gauge (gauge length 2mm) (FLA-2- 
11, Tokyo Sokki, Tokyo, Japan) was bonded to the bottom 
of the midspan. The crack shear displacement and loading 
point displacement were measured by the CSD gauge and 
dial gauge set on the upper cantilever of the specimen and 
below the loading nose, respectively. Two sheets of Teflon 
film of 0.5 mm thickness were inserted between the crack 
surfaces to reduce the friction between the upper and lower 
cantilever beams. This specimen was supported with the 
span whose lengths are shown in Table 1, and the load was 
applied to the midspan of the longitudinal-radial plane not 
to exceed the strain rate of 0.015/min at the midspan. The 
loading speed corresponding to the span is shown in Table 
1. The load R, crack shear displacement d. and longitudinal 
strain at the bottom of the midspan s x were measured si¬ 
multaneously and recorded in a data log at intervals of 2 s. 
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Table 1. Crosshead speed corresponding to each span and depth 


Span 2 L (mm) 

Crosshead speed (mm/min) 

150 

0.3 

200 

0.5 

250 

1.0 

300 

1.0 

350 

2.0 

400 

2.0 

450 

3.0 


The critical load, which is regarded as the load at initiation 
of crack propagation, was determined to be the load be¬ 
tween the straight-line segments through the prelinear and 
postlinear portions of the load-CSD curves . 6-7 Equation 
( 8 ) suggests the difficulty of determining the critical load 
from the load-longitudinal strain relation because this 
relation has no crack length and is not sensitive to crack 
propagation. The load-strain relation was used to calculate 
the energy release rate; hence, the critical load was 
measured from the load-loading period relation. Prior to 
the tests, straight lines were drawn in the crack-free region 
of an LT plane at intervals of approximately 0.03L per¬ 
pendicular to the crack for the specimen. A load was 
applied until the splitting of a marker line was observed. 
This process was performed under a magnifying glass. The 
specimen was then unloaded until a certain load was 
attained. The specimen was reloaded until the splitting of 
another marker. This load-unload process was repeated, 
and the fracture toughness at the beginning of crack propa¬ 
gation G llc and that during the crack propagation G IIR were 
calculated from Eqs. (7) and (10). As mentioned above, the 
crack length was evaluated by a (obtained by observation) 
and by a w and a c [calculated using Eqs. (13) and (17), re¬ 
spectively]. Plural /^-curves corresponding to each span/ 
depth ratio 2Ll2h were obtained from the various fracture 
toughness-crack length relations and were compared with 
each other. 


Compression tests 

To evaluate the crack length by Williams’s end correction 
theory, a w . Young’s moduli E x and E y , and the shear modu¬ 
lus G xy were determined by compression tests. Short-col¬ 
umn specimens whose dimensions were 60 X 30 X 30 mm 
were prepared. When measuring E x and E y , the long axis of 
the specimen was male to coincide with the longitudinal and 
tangential directions of wood, respectively, whereas when 
measuring G xy the long axis was made to coincide, with its 
direction inclined at 45° with respect to the grain. Strain- 
gauges were bonded at the centers of LT planes, and a 
compression load was applied at a crosshead speed of 
0.5mm/min. Young’s moduli were obtained from the stress- 
strain relation. The shear modulus was determined from the 
following equation. 

° xy = 2(1 + 5 y 45 ) ( 18 ) 
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Fig. 3. Fracture toughness at the beginning of crack propagation G llc 
corresponding to the initial crack length/half-span 2L/2/i. Circles 
and horizontal bars represent the means and standard deviations, 
respectively 


where £45 and v 45 are Young’s modulus and Poisson’s ratio 
in the direction inclined at 45° with respect to the grain, 
respectively. 


Results and discussion 

Figure 3 shows the fracture toughness at the beginning of 
crack propagation G nc corresponding to the span/depth 
ratio 2LI2h. When the span/depth ratio was less than 20, 
the G nc value obtained from the load (P)-CSD (3) relation 
increased with decreasing span/depth ratio. This tendency 
was due to the additional crack shear displacement caused 
by the shearing force . 2 In contrast, G„ c obtained from the 
load (E)-longitudinal strain (ej relation was rather stable 
over all span/depth ratios because the longitudinal strain 
was free from the influence of shearing force. When the 
fracture toughness was obtained from the P-d relation, the 
crack length implicitly contained in Eq. (7) was estimated to 
be longer than the real one because of the end correction 
effect shown as Eq. (11), which indicates that the CSD 
compliance is seriously influenced by this effect. In contrast, 
the real crack length, in which the end correction effect is 
ignored, was used when the fracture toughness was calcu¬ 
lated from the P-e x relation, as in Eq. (10). Because of the 
corrected crack length, the value of G llc obtained from the 
P-d relation was always larger than that obtained from 
the P-s x relation. This tendency was commonly found 
in the fracture toughness during crack propagation G IIR . 
When using the P-e x relation for obtaining the fracture 
toughnesses G llc and G IIR , there is an advantage in reducing 
the influence of shearing force. As shown in Eq. (10), how¬ 
ever, there is a drawback that the end correction effect 
cannot be taken into account from the P-e x relation. Addi- 
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Fig. 4. Crack lengths per half¬ 
span calculated from Williams's 
correction theory and the combi¬ 
nation of load-CSD and load- 
strain compliances corresponding 
to the observed crack length per 
half-span. Open and solid circles 
were obtained from Williams's 
correction theory and the combi¬ 
nation of compliances, respec¬ 
tively; dashed line shows the 
condition that the calculated 
crack length is equal to the ob¬ 
served one 
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tionally, there are two drawbacks when obtaining the value 
°f f'llli from the P-e x relation. The first drawback is due to 
the strain-gauge compliance C s , which is not essentially 
changed by the crack propagation, as in Eq. (9). Insensitiv¬ 
ity to crack propagation causes incorrect estimation of frac¬ 
ture toughness. The other drawback is due to that the crack 
length a should be given by the observation when obtaining 
G iir from Eq. (10). As mentioned, the crack propagates 
obscurely in the ENF test, and it is difficult to determine the 
crack length properly. When determining the fracture 
toughnesses G llc and G IIR , use of the P-d relation by the 
specimen with an appropriate span/depth ratio is recom¬ 
mended. 

Using compression tests, the elastic constants E„ E y , and 
G xy obtained were 12.7, 0.88, and 0.72 GPa, respectively. 
Substituting these constants into Eqs. (12) and (13), crack 
length a w was evaluated by Williams’s end correction 
theory. Figure 4 is an example of the normalized crack 
lengths aJL and aJL calculated with Eqs. (13) and (17), 
respectively, corresponding to the observed crack length 
normalized by the half span a/L. As shown in Fig. 4b, the 
tendencies of aJL-a/L and aJL-a/L relations were similar 
when the span was larger than 350 mm (span/depth ratio 
23.3). This phenomenon indicates that the crack length cal¬ 
culated by the combination of load-CSD and load-strain 
compliances was as effective as Williams’s end correction 
theory. As shown in Fig. 4a, however, the value of aJL did 
not increase by crack propagation despite the increased aj 
L when the span/depth ratio was small. This tendency was 
marked when the span was smaller than 300 mm (span/ 
depth ratio 20). The anomalous behavior of aJL was caused 
by the load-strain compliance C s . Figure 5 shows the nor¬ 
malized strain compliance C s IILh corresponding to the nor¬ 
malized crack length, a/L. Equation (9) suggests that the 
strain compliance C s is independent of the crack length, and 
this independence was verified when the span was larger 
than 350 mm. When the span was smaller than 300 mm, 
however, the value of C s was often increased by crack 
propagation, as in Fig. 5. When obtaining the //-curve by the 
combination of load-CSD and load-strain compliances, an 
appropriate span/depth ratio range exists. 
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Fig. 5. Variance of normalized strain compliance by crack propagation 


Figure 6 shows the //-curves, which are the G IIR -a, a c , and 
a w relations and the G IIR -/da, Zla c , and zla w relations. The 
value of G iir was increased by crack propagation due to 
“fiber bridgings”, which span the upper and lower crack 
surfaces; and the //-curves showed the convexity whether 
the crack length was defined by a, a c , or a w because the 
resistance against the crack propagation is increased by 
the bridgings. 5 As in Fig. 6a, the //-curves obtained using the 
observed and calculated crack lengths were different from 
each other. When the fracture toughness and crack propa¬ 
gation length were used as in Fig. 6b, however, the //-curves 
overlapped. 

Considering the difficulty of observing crack propaga¬ 
tion and the inconvenience of Williams's end correction 
theory, the combined compliance method proposed here 
is appropriate for determining the //-curve. As mentioned 
above, however, the span/depth ratio should be large 
enough to reduce the influence of shearing force and to 
increase strain compliance during crack propagation. 
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Fig. 6. Comparisons of R-curves 
obtained by various crack length 
estimations, a Fracture toughness- 
crack length relation, b Fracture 
toughness-crack propagation 
length relation. Fracture tough¬ 
ness was derived by the load-CSD 
relation. Span was 450 mm 




Aa, Aa c , and 4« w (mm) 


Conclusions 

Using mode II fracture in a spruce specimen, the 
toughness-crack propagation length relation (/(-curve) was 
obtained by the ENF test under the constant loading point 
displacement condition. The following results were ob¬ 
tained. (1) The fracture toughness obtained by the CSD 
compliance was influenced by the span/depth ratio because 
of the shearing force that is inevitable in beam deformation. 
The fracture toughness obtained by the load-longitudinal 
strain compliance was free from the influence of shearing 
force, but it was always smaller than that obtained by the 
load-CSD compliance because of ignorance of the end 
correction effect. (2) The crack length, calculated by the 
combination of load-CSD and load-strain compliances, was 
close to that derived by Williams’s end correction theory 
when the specimen had an appropriate span/depth ratio. (3) 
Using the fracture toughness (obtained by the load-CSD 
compliance) and the crack propagation length (calculated 
by the combination of load-CSD and load-strain compli¬ 
ances), the R-curve can be determined properly under the 
constant loading point displacement condition. 
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